Transport-driven toroidal rotation in the tokamak edge 
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The interaction of passing-ion drift orbits with spatially-inhomogeneous but purely diffusive radial 
transport is demonstrated to cause spontaneous toroidal spin-up to experimentally-relevant values in 
the tokamak edge. Physically, major-radial orbit shifts cause orbit-averaged diffusivities to depend 
on v\\ , including its sign, leading to residual stress. The resulting intrinsic rotation scales with Ti/ Be, 
resembling typical experimental scalings. Additionally, an inboard (outboard) X-point is expected 
to enhance co- (counter-) current rotation. 



Rotation patterns strongly affected by turbulent mo- 
mentum transport are broadly observed in nature, for 
example in atmospheric flows, stellar interiors, and accre- 
tion disks [l| . Laboratory tokamak plasmas are observed 
to rotate toroidally in the absence of applied torque, 
with edge rotation directed with the plasma current (co- 
current), often proportional to plasma stored energy W 
over plasma current /„, and reaching tenths of the ion 
thermal speed v t i @, y| ■ Such intrinsic rotation is of 
practical as well as fundamental interest, since it stabi- 
lizes certain instabilities 0] and contributes to a sheared 
radial electric field E r , believed to suppress turbulent 
transport [B| . fntrinsic rotation is of special importance 
for the next-generation tokamak ITER, since a-hcating 
(nuclear fusion) applies no torque [f3]. 

The intriguing experimental findings have triggered a 
broad theoretical search for the spontaneous rotation's 
physical origins. Although neoclassical (collisional trans- 
port) effects have been considered 0, HI, extensive ex- 
perimental evidence indicates that turbulence dominates 
momentum transport 0, • Numerical efforts have in- 
vestigated turbulent momentum transport in the core, 
both linearly 13, U | and nonlinearly (T2] ■ Models for in- 
trinsic rotation, also primarily core-focused, have treated 
quasilinear approximations (lol . [TH effects of in- 

homogeneity of the confining magnetic field B [ll , 14 1 , 
nonresonant correlations between the fluctuating radial 
E x B drift ve,t and parallel velocity v\\ [l7j], and Stringer 
spin-up type effects 18(. Some scrape-off-layer (SOL) ef- 
fects have been entertained , but without system- 
atic consideration of the confined plasma's response. 

Momentum transport in the tokamak edge presents 
particular challenges for theory. The turbulence is 
strong, with statistics very different from quasilinear es- 
timates [2l[. It is also strongly anisotropic, with parallel 
fluctuation length Lu two orders of magnitude larger than 
the radial length scale L± characterizing toroidal velocity 
and equilibrium plasma variation 2^-24|. Since parallel 
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fluctuation gradients fcii ~ 1/^|| and the corresponding 
forces are accordingly weak, turbulently-accelerated W|| 
and the resulting nondiffusive effects [15l4l7| are smaller 
than simple diffusive momentum transport by k\\L± <C 1 
for realistic edge parameters. Most of these effects are 
further reduced, actually proportional to a "symmetry- 



breaking" (fcn) < (kp 1 ' 2 ~ Hal 

on the scale length of the major radius Rq, the resulting 
momentum transport effects scale relative to simple dif- 
fusion as L±/Rq <C 1 in the edge [H, 14 1. Further, the in- 
teraction of edge and SOL makes the problem inherently 
radially nonlocal. For example, the amplitude of the (un- 
normalized) turbulent fluctuating potential decreases in 
the radial direction on a short length scale ~ L± 
j26| . Given these experimental facts, the present work 
analyzes a simplified, purely-diffusive kinetic transport 
model, setting parallel acceleration identically to zero but 
retaining a model edge and SOL, passing-ion drift orbit 
excursions, and spatial variation of the diffusivity, find- 
ing differential transport of co- and counter-current ions 
to cause residual stress and consequent intrinsic rotation 
levels similar to those seen in experiment. 

Analysis begins with a model axisymmetric drift- 
kinetic transport equation for the ions (c. f. [271 ]) 

d t fi + vd y fi - Sv 2 (sin j/) d x f t - D (y) d x (e~ x d x fi) = 0, 

(1) 

in which fi(x, y,v,t) is the ion parallel distribution func- 
tion normalized to pedestal-top ion density over thermal 
speed rii | p t / vu | P t and v is the parallel velocity normal- 
ized to Ufi|pt, positive for co-current motion. The ge- 
ometry is simple slab plus geodesic curvature drift, with 
uniform poloidal B§, toroidal B^, and total Bq magnetic 
field strength. The radial position x, poloidal position 
y, and time t are respectively normalized to L^, the 
minor radius a, and the ion transit time aBo/ Bgvu\ p t- 
The effects of nonaxisymmetric fluctuations are mod- 
eled with an inhomogeneous turbulent diffusivity, nor- 
malized to L^BgVtilpt/aBo and assumed separable, with 
arbitrary poloidal dependence D(y), decaying exponen- 
tially in x. The dimcnsionless parameter 5 = qp^-pt/ L^, 
with q = aB$l RqBq the safety factor and pi the ther- 
mal ion gyroradius, indicates the passing-ion orbit width 
relative to the radial turbulence inhomogeneity. 5 takes 
values around 1/4 for typical ASDEX-Upgrade (AUG) 
H-mode parameters J^. Collisions arc neglected, a rea- 
sonable approximation if pedestal-top ions escape with- 
out experiencing a collision, roughly for Vut ci < 1, with 
Vu the velocity-dependent ion collision rate and r cr the 
pedestal ion stored energy over the ion heat flux. For 
thermal pedestal-top ions, vut ct takes values around 1 
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for typical H-modcs in AUG, JET, and DIII-D, thus su- 
pcrthcrmal pcdcstal-top ions tend to exit the plasma col- 
lisionlessly while subthcrmal ones do not. Boundary con- 
ditions are fi(-ao,y) -> f l0 (v), /,-(oo,j/) -> 0, fc(x < 
0,j/o) = fi(x < 0,y + 2tt), fi(x > 0,y o ,v > 0) = 0, 
and fi(x > 0,?/o + 2n,v < 0) = 0, with y the poloidal 
X-point angle. Eq. |T]) is invariant to a rigid toroidal 
rotation u r j g , normalized to vuB^/Bq. Although the rel- 
evant general theorem [29| does not directly apply to this 
non-Hamiltonian model, Eq. ([T]) does trivially conserve 
a simplified toroidal angular momentum J(v + u r ig)/i cfo, 
as well as a density J fa dv and energy J (1 + v 2 /2)fi dv. 

Eq. |T]) can be approximately analytically solved w-by- 
w in steady state for both large and small -D c ff, results 
agreeing for D c s ss 1. The solution procedure is briefly 
described here, details given in [3(J. Since v appears 
only as a parameter, w-dependent variable transforms can 
greatly simplify Eq. (JXJ) . First, use new spatial variables 

x = x — Sv (cosy — cosj/o) , (2a) 

V = D7j (v) f D (y') e -H™y'-™y°) d y', (2b) 
J ya 

with Dyo(v) = Jy° +27r D(y') exp(-5v(cosy' -cosy ))dy' . 
Physically, x is a drift-surface label and D y o an orbit- 
averaged diffusivity. Next, for v < 0, take y — > 1 — y. 
Finally, transform from iton = e x ^ 2 , obtaining 

dgfi = (D eS /4) (dlfi - u- l d u fi) (3) 

for fi(u,y,v), in which D cS (v) = D y0 (v)/\v\. Bound- 
ary conditions are now /»(0,y) = fio(v), /j(oo,y) -> 0, 
/i(u < 1,0) = /i(u < 1,1), and /t(u > 1,0) = 0. The 
normalized flux of particles with velocity v through any 
closed poloidal contour, T(v), takes the simple form 

1 r 1 

I» =--D v0 {v)u~ l 8 u fidy, (4) 

evaluated at any constant u < 1. 

Eq. © shows the original problem to reduce to a one- 
parameter family of otherwise-identical differential equa- 
tions. Remarkably, the spatially-constant effective diffu- 
sivity D c ff depends not only on the magnitude of v, but 
also on its sign! Physically, as shown in Fig. [IJ this re- 
sults from the fact that co- (counter-) current ions' drift 
orbits are displaced major-radially outwards (inwards). 
For the typical case of turbulent diffusivities larger at the 
outboard, counter-current ions experience larger orbit- 
averaged diffusivities. Preferentially exhausting counter- 
current ions represents a co-current residual stress, al- 
though momentum transport at any given spatial point 
is purely diffusive. 

In solving Eq. (J3|), a Laplace transform approach sim- 
ilar to [3l| yielded the exact Green's function, 



G{u,i,r) = 





Figure 1. (color online). Co- and counter-current passing ion 
drift orbits over turbulence, plotted for a straight-down (a) 
and inboard (b) X-point. Darker shading indicates stronger 
diffusivity. Co- (ctr-) orbits are displaced major-radially out- 
ward (inward) as shown, regardless of the sign of B$ or I v . 

in terms of which the solution may be written as 

fa(u,y) = f t0 e- u2/D ^+[ f l (t0)G(u : £,y)dZ. (6) 

Jo 

A first-order iterative approximation f- , obtained us- 
ing /i(£,0) -> / i oexp(-^ 2 /L» off ) in Eq. ©, was demon- 
strated to yield an approximate normalized flux r//jo|i>| 

1 /2 

with absolute error strictly less than min(0.58Z? o g (1 — 
e -V-Deff) ( 0.75/D^ 2 ), tight bounds for large D eS . For 
small -D c ff, a two-region solution may be used, represent- 
ing fa with a Fourier series for u < 1 (edge) and Laplace 
transforming for u > 1 (SOL), requiring continuity in /, 
and d u fi at u = 1, except possibly at the single point 
u = 1, y = 0. The resulting edge and SOL ODEs possess 
explicit solutions in terms of modified Bessel functions. 
Slightly generalizing [32|], continuity at the LCFS then 
requires the edge solution to satisfy 

2 Jo vv-y 8 j 

at u = 1. The resulting dense matrix for the Fourier 
coefficients has been solved numerically at various -D e ff; 
retaining 10000 modes in y. 

The two approximate solutions, 

r J D cS / (l + ail^ff + a 2 D eS ) , D eS < 1 
f M \v\ ~ |-i£» eff d u \ u=l (jP/ho) dy, D cS >l' 

(8) 

are well-approximated for all D c g by 

3W-i b (' + g^)"""( 1 + ^)' w 

with ai = 0.8224, a 2 = 0.1763, c 2 = 4, c 3 = -4ai, c 4 = 
4ai + e 4 /( 1+Ql+a2 ) - 5 - 9e" 4 ^, c 5 = c 7 = 0, c 6 = 8e" 47 , 



Figure 2. (color online), (a) normalized flux as a function 
of the ^-dependent effective diffusivity, with uniform approx- 
imation (thin solid black), small- D c f{ approximation (thick 
dashed red), large- D c g approximation (thick solid green), and 
large- D e g error bounds (thin dashed green), (b) speed distri- 
bution of particle, momentum and parallel heat fluxes, as- 
suming a Maxwellian at the inner boundary. 



eg = e -47 , and Eulcr's constant 7 f» 0.5772. The second 
approximation is used for the simplified explicit forms in 
Eqs. (fTTj) - (f*P2"l) and corresponding plots, the first for all 
other plots. The results of Eqs. © and (O are plotted 
in Fig.[^a), along with the large- D e g error bounds. 

Explicit forms for normalized particle, momentum, 
and parallel heat fluxes may now be obtained for any 
specified fio(v) and D(y). Assuming a Maxwellian 
fio{v) = e~ v / 2 /(27r) 1 / 2 and simple ballooning form 
D(y) = Dq(1 + d c cosy), thus 



D eS (v) = 2nDoe 5vc ° s y° [I (Sv) 



<kli(6v)]/\v\, (10) 



the relevant flux moments may be reasonably approxi- 
mated for small 8 as 
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in which g p (D ) = ln(l + 27rL> /e 7 p 1/2 ). The integrands 
are plotted in Fig. [2Kb), summed over the sign of v. (All 
plots use representative AUG H-mode values D = 0.033, 
d c = 0.8, y Q = -5tt/8, and S = 0.28.) 

The total momentum flux, incorporating w r i g , is just 
VrigT p + II. Since toroidal rotation damping is very weak 
[33J, a vanishing momentum input implies that momen- 
tum flux must also vanish, resulting in the pedestal-top 
intrinsic rotation rate 
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plotted in Fig. EHa). Alternatively, one may balance an 
applied NBI torque with the outward momentum flux 
resulting from the NBI-driven ion heat flux. For zero 



Figure 3. (color online). Normalized intrinsic rotation veloc- 
ity Vmt (a) and unbalanced NBI fraction II/(r p + Qy) (b), 
plotted as functions of drift orbit width S for several values 
of poloidal X-point angle yo, with numerical integrals (thick 
dashed) and analytical approximations (thin solid). 



pedestal-top toroidal rotation, v r i g = 0, one must set the 
unbalanced NBI fraction / un b = {P§bi — ^nbiV-^nbi to 



unb 



f c i; N bi IT_ 

2/nbi B q VtiUTP + i 



'II 



(13) 



in which /nbi is the fraction of heating by NBI, f c is 
the fraction of heat transported by ions, and i>nbi is the 
beam ion velocity. The ratio n/(r p + Q\\) is plotted in 
Fig. [Hb). Since UNBi/««| P t is typically large, /unb may 
be a significant fraction of unity, as observed by 34 1. 

Several features of this solution deserve comment. 
First, the steady-state results given in Eqs. (fT2"j) - (fllS|) 
are due to a balance between large momentum transport 
terms [Fig. [2Kb)], thus are robust. Relaxation time to the 
edge intrinsic rotation profile should occur roughly on an 
ion transport time through the pedestal, ~ r cr . The v- 
asymmctric diffusion II is independent of the toroidal ve- 
locity and its radial gradient, thus —II represents a resid- 
ual stress. For typical experimental parameters, it acts 
in the co-current direction with experimentally-relevant 
magnitude. The dimensional intrinsic rotation predic- 
tion, given here for small Dq, is 
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The 1/Bg dependence corresponds to experimentally- 
observed 1/Ip scalings 0, while proportionality to 
Ti | p t provides an alternative explanation for recent obser- 
vations [24], Co-current spin- up at the L-H transition 
0,13] is expected due to the increase in Ti| pt and proba- 
ble decrease in L^. The predicted dependence on X-point 
poloidal angle has yet to be experimentally tested. The 
physics presented here may also have implications for in- 
ternal transport barrier (ITB) rotation: for outboard- 
ballooning and radially-increasing diffusivity (as outside 
an ITB), the asymmetric diffusivity results in a counter- 



current core rotation increment, consistent with [361.137 1. 

Simplifications used in this model must be kept in 
mind. The presented calculations omitted both the V-B 
drift and the radial electric field E r , outside the latter's 
contribution to i> r i g . While the VB drift has little effect, 
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a uniform uncanceled poloidal E x B drift of magnitude 
approaching vuBg/Bo can represent a nonncgligible co- 
(counter-) residual stress for outwards (inwards) E r , a 
transport effect due to a shifted relation between D y0 and 
D e g [30l | , possibly connected with observations of a favor- 
able/unfavorable X-point dependence of L-mode toroidal 
rotation [2(|[37j. Treatment of sheared E r effects or elec- 
trostatic confinement [38| would require nontrivial exten- 
sions to the theory, as would retention of collisions. Di- 
rect collisional effects on the rotation-driving flux II may 
often be small, since II results dominantly from ions that 
are slightly superthermal at the pedestal top [Fig. [2fb)], 
thus very superthermal at the LCFS, with an accordingly 
low collision rate. However, lower-energy ions may affect 
both E r and the rotation saturation u r i g r p . Finally, re- 
call that the turbulence parameters are taken as an input 
to the present model, not calculated self-consistently. 

In summary, radial displacements of passing-ion or- 
bits and typical tokamak-edgc turbulence inhomogeneity 
are shown to result in orbit-averaged diffusivitics that 
depend on the sign of v\\. Even in the absence of non- 
diffusive effects, this results in residual stress and cor- 
responding intrinsic rotation at experimentally-relevant 
levels. The rotation is co-current for typical H-mode pa- 
rameters and scales with Ti\ pt /Bg, in agreement with ex- 
perimental observations. 
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